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Abstract. Galaxies in clusters and groups moving through the intracluster or intragroup medium (abbreviated 
ICM for both) are expected to lose at least a part of their interstellar medium (ISM) by the ram pressure they 
experience. We perform high resolution 2D hydrodynamical simulations of face-on ram pressure stripping (RPS) 
of disk galaxies to compile a comprehensive parameter study varying galaxy properties (mass, vertical structure 
of the gas disk) and covering a large range of ICM conditions, reaching from high density environments like in 
cluster centres to low density environments typical for cluster outskirts or groups. We find that the ICM-ISM 
interaction proceeds in three phases: firstly the instantaneous stripping phase, secondly the dynamic intermediate 
phase, thirdly the quasi-stable continuous viscous stripping phase. In the first phase (time scale 20 to 200 Myr) 
the outer part of the gas disk is displaced but only partially unbound. In the second phase (10 times as long as 
the first phase) a part of the displaced gas falls back (about 10% of the initial gas mass) despite the constant 
ICM wind, but most displaced gas is now unbound. In the third phase the galaxy continues to lose gas at a rate 
of about 1 Mq yr~^ by turbulent viscous stripping. We find that the stripping efficiency depends slightly on the 
Mach number of the fiow, however, the main parameter is the ram pressure. The stripping efficiency does not 
depend on the vertical structure and thickness of the gas disk. 

We discuss uncertainties in the classic estimate of the stripping radius of lGunn fc Gotlj lll972l) . which compares 
the ram pressure to the gravitational restoring force. In addition, we adapt the estimate used bv lMori fc Burker3 
Boon) for spherical galaxies, namely the comparison of the central pressure with ram pressure. We find that the 
latter estimate pred icts the radius and m ass of the gas disk remaining at the end of the second phase very well, 
and better than the lGunn fc Got"3 ilQT^T) criterion. 

From our simulations we conclude that gas disks of galaxies in high density environments are heavily truncated 
or even completely stripped, but also the gas disks of galaxies in low density environments are disturbed by the 
flow and back-falling material, so that they should also be pre-processed. 
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1. Introduction 

The differences observed between disk galaxies in clusters 
and in the field indicate that the evolution of galaxies 
is not only ruled by interior processes, but also strongly 
by the galaxy's environment. Some cluster galaxies show 
distortions in their gas components as well in as their 
stellar disks. These examples can be explained by tidal 
interactions with some fellow cluster mem bers (mergers 
or harassment, see review bv lMihoslEooi and references 
therein). However, a number of galaxies have rather nor- 
mal stellar disks but truncat e d gas disks (e.g. N GC 4522: 
iKennev fc KoopmannI Il999l 1200 JVollmer et al. 2004ll 
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NGC 4548: IVollmer et alJll999t NG C 4848: IVollmer et alJ 
l2001a|) . ICavatte et all l)l990lll994|) observed the brightest 
spiral galaxies in the Virgo cluster and found the trend 
that the galaxies closest to the cluster centre have the 
smallest HI disks. In general, cluster spirals are HI de- 
ficient compared to their field relatives. They are more 
deficient near the cluster centre ijSolanes et alJl200lj) . but 
deficient galaxies can be observed ou t to large distances ( 2 
Abell radii) from the cluster centre ijSolanes et alJl200]|) . 
HI deficient galaxies tend to be on more radial orbits which 
lead them deep er into the cluster. iKoopmann fc KennevI 
l|l998l l2004allJ) studied the star formation rates in Virgo 
spirals and found that also the star formation of many 
cluster members is truncated, that means it is suppressed 
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significantly in the outer regions, but normal or even en- 
hanced in the central region of the galaxy. 

The combination of an intact stellar disk but truncated 
gas disk can hardly be achieved by tidal interactions, as 
this usually also disturbs the stellar component. But even 
if a cluster galaxy does not interact with another galaxy, 
it still moves through the intra-cluster medium (ICM) and 
experiences a ram pressure Pram ~ Picm • '^icm; ^^^^ 
pends on the ICM density picM and the relative veloc- 
ity 'ViCM between the ICM and the galaxy. What is ex- 
pect ed to happen t o suc h a galaxy? For a face-on mo- 
tion, [Cun^^Got^ l|l972() suggested to compare Piam and 
the gravitational restoring force per unit area ^Sgas(r), 
which depends on the gradient of the galactic potential 
$ in the direction perpendicular to the disk (z-direction) 
and the surface density of the gas disk Sgas- At radii r 
where 



9$ 

Pram > "Q^^^^ ^gas(^)i 



(1) 



the gas will be pushed out of the disk. Applying this esti- 
mate to typical disk galaxies in clusters, galaxies passing 
near the cluster core should get stripped of large parts of 
their gas disks. In addition to the ram pressure stripping 
(RPS), a remaining inner gas disk should be compressed. 

As the common interest concentrates on the effect on 
the galaxy, numerical simulations of RPS are usually done 
in the rest frame of the galaxy. Here the relative motion 
translates into an ICM wind. 

First attempts to simulate RPS of disk g alax- 
ies have been made by lFa,rouki Sha,nircl ijlflSfll) and 
iTovama fc Ikeucliil l|l98fl() . but these works were limited to 
low resolution and short run times due to weaker compu- 
tational power. Late r studies with bett e r resolution (e.g. 
SPH simulations by lAbadi et all Il999l ISchulz fc Str^ick 
22Jil and 3D Eulerian hydro simulations bv lOuilis et al 



20nf)|) performed a few simulations of massive disk galax- 
ies in constant mostly face-on winds. The simulated winds 
were representative for the conditions in cluster centres, 
the galaxies suffered substantial gas loss. In general, previ- 
ous works found that the gas disks are stripped on a short 
time scale (few lOMyr to few 100 Myr), and that the ra- 
dius of the remaining gas disk is in g ood agreement with 
the estimate derived f r om E q.[Tl Both lOuihs et alJ l)2000|) 
and ISchulz fc Struck! 1I2OO I') used longer run times and 
found a long phase of turbulent viscous stripping following 
the rather short RPS phase. This viscous stripping works 
on galaxies that are oriented inclined or edge-on towards 
the ICM wind as well as on face-on oriented ones. The 
RPS phase seems to take longer for inclined cases, but it 
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are not homogeneous but porous and tried to model this 
feature by holes in the disk. In such cas es the turbu- 
lent vis cous stripping was more effective. IVollmer et al.l 
2001bl) applied a sticky particle code to simulate the pas- 



sage of disk galaxies through clusters, stressing the fact 



that in reality the galaxies experience a rapidly chang- 
ing ICM wind instead of a constant one. Consequently, 
the stripping is a distinct event. Moreover, during this 
event not all displaced gas is unbound, and up to 10% of 
the ii iitial gas mass fall back into the disk. In furthe r pa- 
pers (IVollmer et al.ll20fl(ll2001alll999HVollrneill2003t) this 
group tried to explain the features of individual galax- 
ies by comparing simulation results with observations and 
specifying the ram pressure, inclination angle and t he mo- 
ment of the stripping event. iMarcolini et al.l l(200.'^ simu- 
lated RPS of disky dwarf galaxies in poor groups, where 
the ICM densities and hence ram pressures are weaker 
than in clusters. 

Most previous works performed a few simulation runs 
with massive galaxies in winds representative for cluster 
centres and found that for these cases the gas loss can be 
predicted by Eq. (at least for face-on winds) . Massive 
galaxies in weak winds corresponding to conditions in clus- 
ter outskirts were not simulated, as in general it is assumed 
that weak winds do not have important effects. Here we 
present a comprehensive parameter study of massive to 
medium-mass disk galaxies in face-on ICM winds which 
are representative for the whole range from cluster cen- 
tres to cluster outskirts and groups. We want to check if 
we can confirm the analytical estimate for this large range 
of ram pressures. Moreover, the analytical estimate incor- 
porates only Pram and does not distinguish with which 
combination of picM and ?;icm this Pram was achieved. 
We want to test systematically if the analytical estimate 
misses some dependence on Mach number, especially as 
viCM is expected to be in the transonic range. As a second 
point, we want to investigate the dependence of RPS on 
galactic parameters. Concerning the dependence on over- 
all galaxy mass, we perform simulations for a massive and 
a medium size disk galaxy. In addition, in this paper we 
pay special attention to the vertical structure and thick- 
ness of the gas disk, as the analytical estimate (Eq.QJ does 
not distinguish such cases, either. The influence of an in- 
homogeneous gas disk will be studied elsewhere (Roediger 
& Hensler, in preparation). 

The outline of this paper is as follows: In Sect.|21we 
briefly introduce our code. In Sect. 13 we describe our 
initial model. Section ^ gives an overview of the simula- 
tions we performed and explains some analysis techniques. 
Section 13 presents the results and in Sect. we discuss 
them and draw conclusions. 

Throughout the paper we will use r and z for the radial 
and axial cylindrical coordinates and R for the spherical 
radius. 

2. Numerics 

To answer our questions we perform simulations using a 
2D Eulerian hydrodynamics code (assuming cylindrical 
symmetry ) deve loped from an original implementation of 
iRozvczkal l|l985|) . The restriction to 2D reduces the com- 
putational effort significantly, allowing to perform a large 
parameter study. The code is very similar to the one de- 
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scribed inlNorman fc \yinkleil l)l986|) and the ZEUS code 
ijStone fc NormarJll992|) . Detailed descriptions of the nu- 
merical schemes can be found in these papers. A special- 
ity of these codes is the staggered grid: scalar quantities 
(like density and pressure) are defined at the cell centres, 
whereas the components of vector quantities (like veloc- 
ity and momentum) are defined at the cell walls. For the 
computat ion of the gravitational po tential an ADI solver is 
used (see lNorman fc Winkler! llQS^) . The code uses cylin- 
drical coordinates {z, r) with the symmetry axis at r = 
and the inflow boundary at z = 0. 

The natural viscosity is neglected. However, the code 
makes use of an artificial viscosity in order to handle 
shocks correctly. Without that jump conditions and the 
propagation velocity of shocks could not be reproduced 
correctly. The artificial viscosity is designed to be non- 
zero in compressed regions only. It spreads shock fronts 
over a few cells. Only the diagonal elements of the viscos- 
ity tensor Q are non-zero. Due to the cylindrical symme- 
try also the component Q*^*^ is zero. The two remaining 
components take the form 

Q^^(z,r) = p{z,r) ■ A,viz,r) ■ 

i-CiCs + C2 ■ min{A,i;(z,r),0}], (2) 
Q"^{z,r) = p(z,r) • A^f(z,r) • 

[-Cic + C2 ■ min{A,i;(z, r), 0}], (3) 

where Cs is the adiabatic sound speed, p{z,r) the local 
mass density, Azv{z^r) and Arv{z,r) are the local gradi- 
ents of the velocity field in z- and r-direction, respectively. 
T he constant C^ sets a linea r term, which is, according 
to 'Nor man fc WinkW ('1986'), used rarely "and then only 
sparingly to damp oscillations in stagnant regions" . For 
the application here it is also set to zero. The constant C'2 
sets the quadratic term, it gives roughly the number of 
cells over which a shock front is spread. It is usually set to 
values between 1 and 4. Here we use C2 = 2 in agreement 
to previous work and in agreement with results from nu- 
merical tests. A concern may be that due to the artificial 
viscosity the code cannot handle hydrodynamical insta- 
bilities like the Kelvin-Helmholtz instability correctly. We 
have tested the influence of the viscosity in our simula- 
tions by varying the parameter C2 (see Appendix^ and 
found that the differences are negligible. 

We do not yet include radiative cooling, because to 
do so properly we would need to resolve the multiphase 
characteristic of the interstellar medium (ISM), either by 
a multiphase code or by true spatial resolution. However, 
this is not our aim here. In previous works that included 
radiative cooling a lower limit for the ISM temperature 
was used to prevent cooling to temperatures below 10^ or 
10^ K, which prohibits to study the interesting tempera- 
ture range of th e warm HI gas and also molecular clouds. 
iMarcolini et al.l |2003) tested in how far including radia- 
tive cooling leads to different results concerning the gas 
loss from the galaxy and found no difference whether they 
include cooling or not. 



The 2D treat ment restricts the simu lations to face-on 
winds, however, IMarcolini et al.l l)2003j) showed by com- 
parisons with 3D simulations that the 2D simulations can 
handle face-on winds correctly. What we can learn for 
other inclinations from the pure face-on ones will be dis- 
cussed in Sect. 

In contrast to most applications of ZEUS-like codes 
for studies of a flow past an object, we do not choose 
an open boundary condition for the boundary at Tmax- 
Actually the conditions termed "open boundary" only 
mean a semi-permeable boundary, where the outgoing flow 
is extrapolated over the boundary, but no flow can enter. 
For consistency with the generally established expressions 
we continue to use the term "open boundary" for the con- 
ditions just explained. We found that for long run-times an 
open boundary at Tmax leads to slow ongoing loss of ma- 
terial through that boundary, which gradually decreases 
the density of the wind. As this behaviour is not desirable, 
we implemented a "solid wall" at rmax, which means that 
the velocity perpendicular to the boundary is set to zero, 
nothing can cross the wall. The boundary at z^ax is open 
in the sense explained above. 

3. Initial Model 

3.1. ICM 

In the beginning of the simulations the galaxy is embedded 
in a homogeneous ICM at rest. ICM densities and temper- 
atures cover the range from the centres to the outskirts 
of clusters (10^^ cm^"^ to 10^^ cm^"^ for particle densities 
and 10^ K to a few times 10^ K for temperatures). For 
further details see Sect. 14.31 

3.2. Galaxy model 

We use two galaxy models, one resembling a massive 
galaxy with a rotation velocity at the level of 200 kms~^, 
the other resembling a medium mass galaxy with a ro- 
tation velocity of about 150 kms~^. Both galaxies are 
composed of a dark matter (DM) halo, a stellar disk, a 
stellar bulge and a gas disk. All components but the gas 
disk are static, they provide the gravitational potential of 
the galaxy. The specific model used for each component 
is introduced below, the parameters for all components 
(masses and scale lengths) are summarised in Table ^ 
Fig.^shows the decomposition of the rotation curves into 
the contributions of the components for both the massive 
and the medium mass galaxy. Fig. El demonstrates the cu- 
mulative masses for the stellar disks and DM halos for 
both cases. 

3.2.1. Stellar disk 

The density of the stellar disk follows an exponential 
model: 
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Fig. 1. Rotation curves and contributions of the single components for the massive (a) and the medium size galaxy 
(b). 



Table 1. Model parameters for the massive (standard) 
and the medium mass galaxy. See also text for further 
explanations. 
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The value given is b^^sir — 11.5kpc). For the simple expo- 
nential disk this value is valid for all r; for the flared disk 
b{r) varies according to Eq. ^| 
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Fig. 2. Cumulative masses (mass included in a sphere of 
radius i?) for the stellar disk and the DM halo for the 
massive and the medium size galaxy. 



of the stellar disk). The potential needs to be computed 
numerically. 



3.2.2. Stellar bulge 

The bulg e is modelled following the suggestion of 
iHernauisq l)l993|) . We assume a spherical bulge, which is 
described by 



Pbulgc(^) 



-^^bulgc -Rbulgc 

27r R{R + i?buigc) ' 



(5) 



where i?buigc is the scale radius of the density distribution 
and Mbuigc the total mass. 

3.2.3. Dark matter halo 



For the DM halo the model bv lBurkerd l|l995(l is used. The 
potential for this density dist ribution can be calcula ted 
analytically, it can be found in lMori &: Burker3 l)2000|) . 



p{R) 



DM 



(R + Rum){R^ + Rum) ' 



(6) 



First this model was derived from rotation curves of dwarf 
disk galaxies, but later on this form of DM halo could also 
be confirmed for normal disk ga laxies l|Salucci fc Burker3 
I2OOOI: iBorriello fc Saluccil lioOljl . According to observa- 
tions, the core radius Rum, core mass Mq (mass inside 
^dm) and central density puM follow the scaling relations 



Pdm = 3-10 



M 



kpc 



-2/3 



g cm 



(7) 



The parameters are the total mass M, , the scale radius a* 
and the scale height 6* (the asterisks mark the parameters 



so that Burkert-halos are essentially a one-parameter fam- 

iiy- 
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3.2.4. Gas disk 

We use two different models for the gas disk. In general, 
for both the density distribution can be written as 



p{z,r) 



-Mgas 



Do a?] 

and the surface density as 
Afpas f r 



exp 



(8) 



(9) 



where as usual Afgas is the total mass of the gas disk and 
Ogas, ^gas are the scale radius and the scale height, respec- 
tively. The parameter as is the scale radius for the surface 
density. 

For option one we choose 5gas('") = &o = const.. This is 
again a usual exponential disk like in Eq. ^ For this case 
also flgas and as are identical. 

The second option is a flared disk with 

bgas{r) = 6oCxp . (10) 

Here Ogas and as are not identical, but they relate with 
the parameter rt according to: 



1 

n 



1 

as 



OSflgas 

n = — — — ■ 



(11) 



To have an unperturbed flow past the galaxy, there 
must be enough space between the outer edge of the gas 
disk and the r^ax boundary of the computational grid (we 
demand that the disk radius is smaller than about 1 /4 of 
fmax)- Due to the large as the gas disk would extend to 
large r, leaving not enough space between the disk edge 
and the grid boundary, although the grid is already large 
(see Sect. 14. To prevent such problems, we cut the gas 
disk gradually to a finite radius of 26 kpc. If the disk gets 
stripped at radii smaller than this, also the parts outside 
26 kpc would be stripped, so this does not bias the results. 
Due to this cut the true mass of the gas disk Afgas is a bit 
lower than what is set by the parameter A^gas in Eq. |S1 
The values given in Table are the true masses Afgas. 

For stability reasons the gas disk needs to be in pres- 
sure equilibrium with the surrounding ICM and in rota- 
tional equilibrium with the galactic potential. Therefore, 
after setting all densities and calculating the gravitational 
potential of the galaxy, we solve the hydrostatic equation 
in ^-direction. Thus knowing the pressure distribution we 
compute the rotation velocity for the gas disk so that 
in radial direction the gravitational force, the centrifu- 
gal force and the force due to the radial pressure gradient 
are in equilibrium (the pressure force is rather small in 
radial direction) . The resulting rotational velocities along 
the galactic plane for the massive and the medium galaxy 
are shown in Fig. ^ Density contours, the temperature 
distribution and the rotation field for both the normal ex- 
ponential gas disk and the flared one can be seen in Figs. 13 
and0I respectively (both Figs, are for the massive galaxy). 
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Fig. 3. Contour plots for the density (left), temperature 
(middle) and rotation field (right) for the massive galaxy 
with a simple exponential disk. Both the density and the 
temperature contours have a logarithmic spacing. For the 
density plot there is one contour line every half order of 
magnitude; the inmost contour at 10~^'*gcm~^. For the 
temperature the spacing of the contour lines is one or- 
der of magnitude; the contours are labelled in units of 
logio(T/ K). In the velocity field the contours are labelled 
in kms"'^. Please note that the disk appears thicker be- 
cause the axes have different scales. 
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Fig. 4. Same as Fig.Ofor the massive galaxy with a flared 
disk. 



In addition. Fig. [S] provides an overview of profiles of the 
density, surface density, pressure and temperature along 
the galactic plane. Also the rotation curves are repeated 
as they differ slightly for the exponential and the flared 
disk. Fig. |S1 also shows profiles for a thicker and a thin- 
ner (larger and smaller scale height) exponential disk (see 



Sect. ^ 

We want to stress that already the high external pres- 
sure has a strong influence on the gas disk. If the galac- 
tic potential and the density of the gas disk are given 
and the gas disk shall be in pressure equilibrium with the 
surrounding ICM, then the pressure in the disk increases 
with increasing ICM pressure in the way demonstrated in 
Fig. El for the exponential disk. As explained above, for 
the given gas disk density and galactic potential the pres- 
sure in the gas disk was computed assuming hydrostatic 
equilibrium with the ICM in z-direction. The resulting 
pressure and temperature profile inside r = 20 kpc for the 



6 Elke Roediger and Gerhard Hensler: 




- 100 - ; - 

50 - i - 

I 1 1 1 1 lJI 

5 10 15 20 25 30 

r / kpc 

Fig. 5. Profiles along the galactic plane for the massive 
galaxy. The panels show profiles for the density surface 
density S, temperature T and the rotation velocity Vrot 
of the gas disk; all quantities are shown for the case of a 
normal exponential disk (thick solid lines), thick and thin 
exponential disk (thin solid lines) and a flared disk (thick 
dashed lines) (see Sect. 14. 4.1)1 . The surface density is the 
same for all cases. We also show the density profiles for the 
stellar component (disk+bulge) and the DM halo, and the 
surface density of the stellar component. The thick lines 
correspond to the cases shown in Figs. 01 and 0] 



pressure stripping of disk galaxies 
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Fig. 6. Profiles of the density p, pressure p and temper- 
ature T along the galactic plane, as they result from the 
hydrostatic equilibrium for different ICM densities. The 
profiles are shown for the case of an exponential gas disk. 
The density profile for the disk was fixed, it does not 
change with changing ICM density. The ICM densities, 
pressures and temperature can be read from the plot on 
the rhs of the dashed vertical line (the ICM temperature 
of 4.385-10'' K was the same in all cases). This plot demon- 
strates how everywhere except in the very centre the pres- 
sure in the gas disk adapts to the ICM pressure (see text 
for a discussion). 

lowest ICM pressure (solid lines in Fig. O are represen- 
tative for the case of an isolated galaxy (one not embed- 
ded in an ICM). From the pressure profiles in Fig. El it 
is already obvious that e.g. for an ICM with (particle) 
density of about 10""^ cm~"^SlO~^^ gcm~^ and tempera- 
ture TicM ~ 4 • 10^ K, which corresponds to conditions 
in the cluster centre, its pressure is higher than the pres- 
sure in the disk of an isolated galaxy at radii larger lOkpc. 
For this example, in case of hydrostatic equilibrium with 
the ICM, the resulting pressure and temperature profiles 
agree with the profiles in an isolated galaxy only in the 
inner 3kpc. Outside this radius the pressure and tempera- 
ture of the disk are much higher than in an isolated galaxy. 
However, as we want to start with a stable model, we use 
the hydrostatic solutions as initial models. In how far this 
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unrealistic pressure and temperature profiles influence our 
results will be discussed in Sect.|Hl 

4. Simulations 

The main purpose of this paper is to investigate how 
the mass loss from the gas disk due to RPS depends on 
the strength of the ICM wind, and on the galaxy itself. 
Therefore one needs to vary ICM parameters as well as 
galaxy properties. A further critical point we check is how 
the initialisation of the simulation influences the results. 
In this chapter we list the models we have performed and 
explain our choices of parameters. In addition, we explain 
some analysis techniques we used to obtain the results. 

4.1. Grid size and Preparative Tests 

The numerical grid must be large enough so that the 
boundaries do not disturb the flow. After some tests we 
chose a grid size of 100 kpc x 100 kpc and a resolution 
of 153 pc (650^ grid cells). We tested the influence of the 
resolution by running some representative cases with res- 
olutions of 50 pc, 100 pc and 200 pc and found that the 
results are not sensitive to this (see Appendix Ib)| . The 
galaxy is placed at a distance of Zgai — 40 kpc from the 
inflow boundary. We also tested that the exact position of 
the galaxy does not influence our results as long as it is 
not placed too close to the boundaries. We also checked 
that our results are not influenced by the assumption of 
the static gravitational potential by performing a control 
run which took the evolution of the potential due to the 
changing gas distribution into account. Wc found that our 
results are not sensitive to this, because the gas disk con- 
tributes only a minor fraction to the total potential. 

4.2. Wind initialisation 

Initially all gas is at rest. To start the simulation we in- 
crease the ICM velocity at the inflow boundary and the 
outflow boundary linearly over the time interval tswitch 
from zero to the flnal value wicm- Afterwards the inflow 
conditions are kept at the upstream side boundary, and 
the downstream side boundary switches to the usual open 
boundary condition (if ficM is supersonic, the outflow 
velocity at the downstream side boundary was only in- 
creased up to Mach 0.9 and then this boundary switched 
to the open condition). So tswitch is the time during which 
the simulation is "switched on" . Two contradicting de- 
mands make it difficult to choose this parameter. On be- 
half of a smooth flow ^switch should be as long as possible. 
On the other hand, we expect that the galaxy responds 
to the ram pressure on short time scales (a few 10 Myr, 
see Sect. 15. lH . Therefore tswitch should be rather short 
in order to resolve this time scale. As this choice is not 
straightforward, we run some representative cases with 
^switch = 20, 70, 150 Myr. The shortest iswitch is chosen 
to resolve the expected stripping time scale fScct. 1^1.1(1 . 
The longest tswitch was chosen to be as long as possible 
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Table 2. Parameters for the simulations testing influ- 
ence of tswitch (see Sect. I4.2|) . AU these runs use the mas- 
sive galaxy with an exponential disk with 6gas — 0.4 kpc 
and TicMi, corresponding to a sound speed cicm — 
1000 km s~^ 



Pram 


VlCM 


niCM 




tswitch 


(cm^^km^s^^) 


(kms^i) 


( cm~^) 




(Myr) 


1000 


800 


1.56 ■ 10" 


-3 


20, 70, 150 


1000 


2530 


1.56 ■ 10" 


-4 


20, 70, 150 


100 


800 


1.56 ■ 10" 


-4 


20, 70, 150 


100 


2530 


1.56-10" 


-5 


20, 70, 150 



under the condition that the initialisation must be fin- 
ished before the first stripped material reaches the outflow 
boundary. Stripped material that would reach the outflow 
boundary during tswitch would cause reflection problems, 
as during tswitch the outflow boundary has a fixed veloc- 
ity and is not open, as explained above. The maximum 
tswitch was found by trial and error. The parameters for 
these runs are listed in Tabled 

4.3. Wind parameters 

To investigate the influence of the wind strength we ex- 
pose the massive galaxy with a standard gas disk (ex- 
ponential, 6gas — 0.4 kpc) to different ICM winds. These 
winds resemble conditions reaching from the centres to the 
outskirts of galaxy clusters. 

Typical velocities of cluster galaxies (and hence the 
wind velocities) are of the same order as the sound speed 
of the ICM. The ICM temperatures in clusters are be- 
tween 10'' K and 10^ K. For the standard ICM tempera- 
ture we choose Ticmi = 4.385-10'' K. With a mean molec- 
ular weight of /X = 0.6 for the ionised ICM this temper- 
ature corresponds to a sound speed of 1000 kms^^. For 
some simulations a lower ICM temperature of T'icm2 = 
7icMi/'\/T0 = 1.39 - 10^ K corresponding to a sound speed 
of 1000 kms-Vv^ = 562 kms"^ is used. If not stated 
otherwise, we use Ticmi- ICM particle densities nicM 
range from 10"^ to a few 10~'^cm~^. We additionally in- 
clude extreme cases for the density to be able to scan the 
ICM parameter space systematically. 

For the sake of suitable numbers we will use the specific 
ram pressure 

Pram = "-ICMWicM (12) 

in units of cm^'' km^ s^^ throughout the paper. It relates 
to the real ram pressure Pram through 

Pram = M ' amu • Pram 

= 1.002 - 10-14 erg cm'^. (13) 

cm 3 km s ^ 

In Fig. 13 we show isolincs of Pram in t he (wicm,"-icm) 
plane. Our simulations cover a range from 10 to 
10 000 cm-3km^s-2 in p,am (lO^^^ ^ iQ-iOerg cm^^ in 

Pram ) ■ 
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Fig. 7. Models in the wicm-'^icm plane. Contours of con- 
stant specific ram pressure pram are plotted with a spac- 
ing of one order of magnitude. The lines are labelled with 
Pram in units of cm~^ km^ s~^. The open circles mark the 
(^icmj'^icm) pairs used for the simulations with the mas- 
sive galaxy and Tjcm i • The triangles mark the winds used 
for simulations of the medium mass galaxy and the same 
TicM- The crosses mark the winds used for the massive 
galaxy and ricM2 (see also Sect. 14. 



Table 3. Parameters for the simulations testing influence 
of the vertical structure of the gas disk (see Sect. 14.4.11) . 
Common to all these runs are the massive galaxy, Tjcmi, 
Switch = 20Myr. The "•" means the same value as in the 
previous line. 





'''fcM 






gas disk 


6gas/ kpc 


1000 


800 


1.56 ■ 10- 


-3 


exp 


0.2, 0.4, 0.8 










flared 


0.2, 0.4 


1000 


2530 


1.56 ■ 10" 


-4 


exp 


0.2, 0.4, 0.8 










flared 


0.2 


100 


800 


1.56 ■ 10- 


-4 


exp 


0.2, 0.4, 0.8 










flared 


0.2 


100 


2530 


1.56 ■ 10' 


-5 


exp 


0.2, 0.4, 0.8 










flared 


0.2 



The units of pram, vicM and nicM are the same as in 
Table H 

^ The value given is b{r — 5 kpc) . For the simple exponential 
disk this value is valid for all r; for the flared disk 6(r) 
changes according to Eq. 1101 



flared gas disks with different 6gas to some representative 
ICM winds. The simulations done for this point of interest 
are listed in Table 13 



We choose pairs of (wicmj'^icm) so that every Pram is 
covered by four different combinations of wicm and rijcM, 
two in the subsonic and two in the supersonic regime. In 
this way we can see if the success of ram pressure stripping 
depends on Pram alone or if it also depends on the Mach 
number. The pairs chosen for the simulations with the 
massive galaxy and Tjcmi are marked by open circles in 
Fig.d 

To study of the influence of the Mach number in more 
detail, the simulations are repeated with ricM2 (marked 
by crosses in Fig. [T)) . We choose Ticm , vicm and tt-icm so 
that the velocities wicm = 253, 450, 800, 1423, 2530 km s^^ 
correspond to Mach numbers 0.253,0.45,0.8,1.423,2.53 
for TicMi, and to Mach numbers 0.8,1.423,2.53,4.5 for 
TicM2, respectively. In this way many cross-comparisons 
are possible. E.g. the stripping effect of a wind with a ve- 
locity of 800 kms^^ can be studied once for the case this 
velocity corresponds to Mach 0.8 and once to Mach 1.423. 

4.4. Influence of the galactic parameters 

In this paper we focus on two aspects of the galaxy. First 
we address the question whether the stripping efficiency is 
determined solely by the gas surface density as suggested 
by the Gunn & Gott criterion (Eq. or if the vertical 
structure of the gas disk plays a role. The second point of 
interest is the influence of the total mass of the galaxy. 

4.4.1. Influence of the vertical structure of the gas disk 

To study the influence of the shape and thickness of the 
gas disk on the RPS efficiency we expose exponential and 



4.4.2. Influence of galaxy mass 

In order to check our results for a smaller galaxy (see right 
column of Table ^ we rerun the simulations with different 
wind parameters using Ticmi- The wind velocities and 
particle densities used here are marked in Fig. [3 by small 
open triangles. 

4.5. Analysis techniques - gas disk mass and radius 

For all simulations we need to calculate the mass and ra- 
dius of the remaining gas disk. Concerning the fate of the 
stripped gas, with our code we can trace how much of the 
stripped gas falls back to the disk and how much is lost 
permanently. 

4.5.1. Definitions 

We use the following definitions: 

The original disk region is defined as the region where 
P > PiCM at the start of the simulation. 

Galactic gas is gas that has been inside the original disk 
region initially. The method used to trace this gas is 
described in Sect. 14. 

We measure at each time t: 

(a) mass Morig(i) of galactic gas inside the original disk 
region, 

(b) mass Mcy\{t) of galactic gas inside a fixed cylinder 
centred on the galactic plane with r < 26 kpc and 
\z - 2gai| < 5 kpc, 
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(c) mass Mbnd {t) of galactic gas bound to the galaxy po- 
tential. Bound gas is identified by a negative total en- 
ergy density Ctot = Spot + etherm + Crot + Ckin (potential, 
thermal, rotational and kinetic energy density). 

(d) mass Mfb(t) of galactic gas that has left the original 
disk region at some moment but is now back in the 
original disk region - that is the amount of fallen back 
gas (see Sect. I4.5.'2)l . 



The radius of the remaining disk r^ii^]^{t) is determined 
by scanning along the disk plane, starting at r = 0. The 
first cell in which the density drops below 10~^^gcm~^ 
defines the radius of the gas disk (see Fig. ll9|) . This seems 
rather arbitrary, but we tested several other versions, and 
this one turned out to be the most appropriate and rep- 
resentative one. 



4.5.2. Colouring technique - tracing the galactic gas 

To follow the fate of the galactic gas we us e the technique 
of "colou ring" the gas. For its basic version (Severing" '19951 
iLohmann 2000') one has to do the following: Define an ad- 
ditional array COLOURl. At the start of the simulation, 
copy the density of the galactic gas to this array. During 
the simulation advect this array in the same way as the 
usual density. Thus the value of COLOURl in a certain 
cell tells the density of galactic gas at this position. 

We extended this method to determine how much of 
the stripped galactic gas falls back to the original disk re- 
gion. Therefore gas of colour 1 (galactic gas) that leaves 
the original disk region "changes its colour" to colour 2. 
This is done with the help of a second array C0L0UR2, 
which is set to zero initially. Then at each time step the 
procedure is as follows: for each grid cell outside the origi- 
nal disk region add the value of COLOURl to the value of 
C0L0UR2 in the same cell, then set COLOURl to zero in 
this cell. This causes the stripped gas to change its colour. 
Stripped gas with colour 2 that comes back to the disk re- 
gion does not change its colour again, but keeps colour 
2. Hence integrating over the array C0L0UR2 inside the 
original disk region gives the amount of galactic gas that 
has been outside at some time - which is the amount of 
galactic gas that has been stripped and fallen back Mfb. 
Integrating over COLOURl gives the amount of galactic 
gas that has never been outside the original disk region 



M 



stay ; 

and Affb -I- M^tay 



Moris is the total amount of 



galactic gas inside the original disk region, no matter if it 
has been outside this region in the meantime or not. 



5. Results 

In this Sect, we first make some analytical estimates and 
explain what we expect to happen. After that we present 
the results of the simulations and discuss them. 



5.1. Analytical Estimates 

As discussed by previous authors (e.g. iMori fc BiirkertI 
[2000; Marcoli ni et alJl2003tl . the ICM-ISM interaction can 
be divided into two sub-processes: instantaneous stripping 
and continuous stripping. The first process is determined 
by the force balance between the ram pressure and the 
gravitational resistance. Depending on the strength of the 
ram pressure, the outer gas disk down to a radius rgtrip 
should be pushed away in a relatively short time (see 
Sect. 15.1.11) . In con trast, the con tinuous (turbulent vis- 
cous) stripping (see l^lgs^ ^13) is a phenomenon con- 
nected with the boundary between the ICM and the gas 
disk and hence works on the whole surface of the gas disk. 
It has a rather long time scale (see Sect. 15. 1.211 . Due to 
the different time scales of the two subprocesses the over- 
all process should appear to happen in two phases, first 
the instantaneous stripping, followed by the continuous 
stripping phase. 

5.1.1. Instantaneous stripping 

As mentioned in the introduction, ICunn fc CottI l(l972^ 
suggested to compare the ram pressure and the gravita- 
tional restoring force per unit area to estimate which part 
of the gas disk will be pushed out by the ICM wind. We 
consider this estimate a bit closer here. A crucial point is 
that it assumes that the gas disk is thin enough that its 
mass distribution can be represented solely by its surface 
density Sgas(r). Secondly, the gravitational acceleration in 
z-direction in any position is ^(z,r). So, strictly speak- 
ing. 



|/grav(2;,7-)| 



9$ 

— (z,r)Ega,(r) 



(14) 



is the gravitational force per unit area in z-direction that 
works on the gas disk at radius r if this gas disk is shifted 
to position z in the galactic potential. We point this out in 
such detail, because for a given r, /grav(2, changes with 
z. E.g. in the most simple the accelera- 

tion due to the gravitational potential §f = and so is 
the restoring force. In Fig.|Slwe plot profiles of /giav(-z, r) 
perpendicular to the galactic plane (in z-direction), on 
the downstream side of the galaxy, for different radii. As 
expected, for small radii the restoring force is strongest. 
Interestingly, with increasing r the maximum of /grav(-z) 
is located at increasingly large distances from the galactic 
plane. So for radii > lOkpc the maximum of /grav(-2) is 
outside the initial gas disk (this feature was already men- 
tioned briefly bv ISchulz & Struckll2001^ . Which /g,.av(^) 
should be used for the estimation of rgtrip? Two plausible 
choices are possible: either to use the maximal fgraviz) in- 
side the original gas disk, or to use the maximal /grav(-z) 
at any z. For the comparison of Pram and /grav('') tbe 
two version can lead to different results concerning radii 
^ lOkpc. 

At the stripping radius fstrip, Pram and /grav(?') are 
equal. For smaller r the gas is retained in the galaxy as 
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Fig. 8. Profiles of the gravitational restoring force per unit 
area in z-direction fgra.v{z,r) behind the disk, at different 
radii r. Each profile is labelled with its cutting radius in 
kpc. The label is placed at the maximum of the particular 
profile. Please note that for larger r the strongest gravi- 
tational resistance is far behind the disk. For comparison, 



the right y-axis is labelled in units of cm km^ 

of Pram)- 



(unit 



■grav,max: i" Original 

everywhere 
pressures: exp 
flare 
fit 





Pram=niCM-vfcM i" cm-3(km/s)^ 



Fig. 9. Analytical estimate of stripping radius Tstrip and 
mass A/strip of the remaining gas disk after the instanta- 
neous stripping as a function of pram. The thin dash-dotted 
and short-dashed lines are derived from comparing the 
gravitational restoring force /grav with the ram pressure, 
the dash-dotted line uses the maximum of /grav(-2) inside 
the original gas disk, the short-dashed line is the estimate 
using /grav (2) at any z. The solid and long-dashcd thick 
lines are derived from comparing the thermal pressure in 
the galactic plane with the ram pressure, for the flared 
and the exponential disk, respectively. See also Sect. 15.11 
The dotted line is the fit to the result of the pressure com- 
parison for the exponential disk according to Eq. |^ 



Pram < /grav(?'), for larger r the gas will be lost because 

Pram > /grav(?')- Thc resulting Tstrip aS a function of Pram 

is shown in the left panel of Fig.|51for the massive galaxy 
for both versions of /grav (2)- The right panel shows the 
mass of the gas disk remaining after thc instantaneous 
stripping Mstrip, i.e. all gas at radii smaller than rgtrip- 
The estimates according to the two versions differ only 
for lower Pram, and the difference is small. 

We note that for gas disks with larger scale heights 
the assumption that the mass distribution is represented 



well enough by the surface density alone may be wrong. 
In this case the gravitational restoring force per unit area 
on a gas package of volume AF is 



I /grav (2,7-) I = 



dz^'' AA 
(9$ 

-^{z,r)pg^s{r)A. 



(15) 



where AA thc surface of the gas package perpendicular to 
the wind direction and Az is the length of this package 
in wind direction. To obtain an upper limit, for /Ogas(?') 
the density in the galactic plane is used. It is not obvious 
which value should be used for Az, so that this is rather a 
free parameter (which should be somewhere between the 
grid resolution and the overall thickness of the disk). If the 
scale height 6gas is the same for all radii r, the expression 
in Eq.^lis proportional to the one in Eq.^l but smaller; 
Eg. 1 141 would overestimate the restoring force. With /grav 
from Eq.^1 the predicted rstrip(Pram) has the same shape 
as before (using Ea. ll4|l . but is shifted in horizontal direc- 
tion towards lower Pram- If however the gas disk is flared, 
the expressions in Eos. 1141 andll5larc not proportional and 
the analytical estimates of ?'strip(Pram) would have differ- 
ent shapes. 

iMori fc BiirkertI l(200f]l) suggest an alternative criterion 
for the stripping efficiency for spherical galaxies, namely 
a spherical galaxy will be stripped completely if the ram 
pressure exceeds the thermal pressure in the galactic cen- 
tre. We modify this estimate for the application to disk 
galaxies. Given the radial thermal pressure profile po{r) 
in the galactic plane of an isolated galaxy (as shown in 
Fig. El , the gas will be stripped from radii r where 



Pram > Po(»')- 



(16) 



Again, rgtrip is the radius where Pram = p('')- The result 
from this estimate is also shown in Fig. for both the 
exponential and the flared disk. As the pressure profiles 
for both cases are similar, also the functions rstrip(Pram) 
are similar. The difference between the pressure profiles 
for the different disk scale heights (thick and thin disks) 
is not larger than the difference between the pressure pro- 
files for the flared and the exponential disk (see Fig. . 
Hence the profile for the medium scale height is represen- 
tative for all thicknesses. Interestingly, the shape of the 
function rstrip(Pram) is very similar to the shape of the 
estimate from the Gunn & Gott criterion, but the result 
from the pressure comparison is offset towards smaller ram 
pressures, indicating thatpo(?') < /gravmax('')- This can be 
understood analytically. The pressure in the galactic plane 
at a certain radius Po{r) derived from hydrostatics (for an 
isolated galaxy, i.e. p(±oo) = 0) is 



9$ 

Poir) = I - — {z)p{z)dz. 



(17) 
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The calculation of |/grav| uses the steepest potential gra- 
dient 1 4^ (r) I and the surface density 

\ oz ^ ' \ max 



E(r) = y p{z)Az = 2 J p{z)dz. 

— CO — oo 

With that we can write | /grav ('') | as 



(18) 



l/grav(?-)| 



5$ 

dz 



r) 



p{z)Az. 



(19) 



As in Eas.1171 andll9lall functions are positive in the inte- 
gral range and as I ^ I > 

° ° \ oz I max — 

that 



|/grav(»-)| > 2po(T')- 



|j in this range, it follows 
(20) 



For our galaxy model we have /grav(?') ~ 4po('')- 

The time scale for the instantaneous stripping tgtrip 
can be estimated from the simple uniformly accelerated 
movement. With a constant acceleration a ~ Pram/Sgas a 
body that starts with zero velocity moves the distance s 
during the time 



^strip 



= 56 




(21) 



5kpcy ylO^-^g cm ^ 

\ -1/2 

Pram 

1000 cm-3 km^ s-^ 



1/2 



Myr. 



(22) 



Let us assume that a gas package is stripped off when it 
has reached a distance of 5kpc to the disk plane. For the 
medium galaxy the surface density is about a factor of 0.38 
lower, hence the corresponding stripping time scale for the 
medium galaxy should be about a factor of VO.SS ~ 0.6 
shorter. This estimate is less accurate for small ram pres- 
sures because in that case also the gravitational restoring 
force needs to be taken into account (for small Pram this 
estimate gives a lower limit). For large ram pressures the 
gravitational deceleration can be neglected because it is 
small compared to the acceleration due to Pram- 

5.1.2. Continuous stripping 

According to lNulsenI l)l982|) the maximal mass loss rate of 
a spherical cloud of radius R in an ICM wind with particle 
density nicM and wind velocity wicm is 



M « 20 



R 



20kpc 



) V 10-3 cm-3/ I, 



1000 kms- 







Although this estimate was derived for spheres, also a 
rough estimate for disks can be deduced. Apart from 
the geometry, for the face-on wind flowing past a disk 
galaxy, the radius of the remaining gas disk depends on the 
strength of the wind, as the gas disk is truncated quickly 




100 1000 10000 

Pram=niCM-vfcM cm-^(km/s)^ 

Fig. 10. Analytical estimate of the mass loss rate for the 
continuous stripping. The solid line shows the result if a 
constant radius for the gas disk is assumed, the dashed 
line shows the reduced mass loss rate if rgtrip as a function 
of Pram IS uscd (Eq. 1241) . Both functions are shown for 
vicM = 1000 kms-^. One has to keep in mind that this 
estimate is only rough, as it is derived for a spherical and 
not a disk galaxy. 

by the instantaneous stripping. Over a wide range of ram 
pressures the stripping radius can be fitted by 



rip 



= 1.037- 0.253 log 



10 



Pr. 



cm-3 km^ s-^ 



rip 



for the initial disk radius rp = 24.54 kpc or 

Pram 



kpc 



25.44- 6.2 log 



10 



' km 



(24) 



as is shown in Fig. We can rearrange Eq. ESI to 
M « 20 f —4^ \ X 



R 
20 kpc 

Pram 



103 cm-3 ]jQ-^2 g-2 j V 103 km s 



-1 M, 







yr 



(25) 



and replace R by 7'strip(Pram), thus deriving a rough esti- 
mate for the mass loss rate for our case. Fig. 1101 compares 
the predicted mass loss rates of the continuous stripping 
as a function of Pram for wicm — 1000 km s-^ for a constant 
R on the one hand, and for R — r strip = ''strip (Pram) on the 
other hand. If we take into account that the radius is trun- 
cated by the instantaneous stripping, the mass loss rate 
is reduced significantly. However, this estimate can only 
be rough as it was derived for spherical bodies instead of 
disks. Beyond that the dependence on t;icM (see Eq. EHJ 
introduces a further factor. So we can only estimate that 
the mass loss rate due to the continuous stripping is of the 
order of one Mq yr-^. With this number the time scale 
for the continuous stripping 

Afgas(i = 0) 



M 



(26) 



(23) 



is of the order of several Gyr. 

5.2. General behaviour of the simulations 

To give an impression of the temporal evolution of the 
stripping process in our simulations, we show the gas den- 
sity distribution at several time steps for a strong sub- 
sonic, a strong supersonic and a weak subsonic wind in 
Figs. llllO andlT^ respectively. 
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Fig. 11. Evolution of gas density distribution for 
the massive galaxy in a wind with Tjcmi, Pram = 
1000cni-3km^s-2, Mach 0.8, ^switch = 20Myr. The den- 
sity distribution is shown by contour lines with the inmost 
contour at lO^^'^gcm^^; the spacing of the contours is 
half an order of magnitude. Only contours for p > picM 
are shown. In the shaded regions the gas is still bound to 
the galactic potential. The size of the small rectangle in 
the top left corner corresponds to 10 x 10 grid cells. 
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Fig. 12. Like Fig. ^2 but for a wind with pram = 
1000 cm^'^ km^ s~^, Mach 2.53. Again contours for p > 
PiCM are shown. Here picM is lower than in Fig. II II hence 
there are more contour lines in this plot. As this is a su- 
personic wind, also a bow shock appears. 



ration of this phase depends on the ram pressure. For 
high ram pressures it is short (a few 100 Myr), for low 



We find that the ICM-ISM stripping proceeds in three 
phases: 

Instantaneous stripping: First the outer gas disk, down to 
a certain radius, is bent towards the downstream side 
of the galaxy. This phase can be seen in the first panel 
with t > Myr in Figs. ^2 to ^1 The time scale on 
which the bending takes place is determined by the 
ram pressure. It reaches from a few 10 Myr for high 
ram pressures to about 200 Myr for low ram pressures. 
Also the radius down to which the disk is bent de- 
pends on the strength of the ram pressure. Reasonably, 
for higher ram pressures only a small inner part is re- 
tained, and for low ram pressures only the outer edge 
is bent, but the disk is affected even for very small ram 
pressures. In all cases a large part of the matter that 
is pushed out of its original position is still bound. 

Intermediate phase: The first phase is followed by a quite 
dynamical intermediate phase, during which the bent- 
out part of the disk breaks up. A part of the material 
that is pushed out is now stripped completely and a 
part falls back to the original disk region. This back- 
fall always happens in the lee of the galaxy, where the 
displaced gas is protected against the wind. This phase 
in which bound gas can be found b ehind the gas disk 
was also noticed by ISchnlz Struck. (,200Ll . The du- 




10 15 20 25 



5 1 1 5 20 25 -5 



10 15 20 25 



Fig. 13. Like Fig. ^2 but for a weak wind with pram — 
10 cm~^ km^ s~^, Mach 0.8. Concerning the number of 
contour lines, see the comment at Fig. El 
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Fig. 14. Evolution of density p and surface density S pro- 
files along the galactic plane for to the runs shown in 
Figs. ^2 (left) andllHI fright). For the calculation of S we 
integrate over interval (^gai — 5kpc, Zgai + 5kpc) to avoid 
that stripped material at larger distances to the galaxy 
appears in E. 

ram pressures this phase continues for the rest of the 
simulation (up to 2 Gyr) . 
Continuous stripping: Actually the continuous stripping - 
the peeling off of the outer disk layers of the upstream 
side of the gas disk by the Kelvin-Helmholtz (KH) in- 
stability - works on the gas disk from the beginning 
on. However, the effect of this process is quite small, 
so that it is not distinct before the first two phases are 
finished. This status is only reached for stronger ram 
pressures, as for weaker ones the intermediate phase is 
long. Even though the effect of the continuous strip- 
ping is small, it slowly decreases the mass and the ra- 
dius of the remaining gas disk. 

Typical for the supersonic runs is a stagnant region in the 
wake of the galaxy, where gas with low densities, that is 
still bound to the galaxy, tends to linger (see final stages 
in Fig.lHl). 

In all runs bound gas can be found at quite large dis- 
tances (10 to 15 kpc) behind the galaxy during the in- 
stantaneous stripping phase and the intermediate dynamic 
phase. For lower ram pressures some gas packages are still 
bound at distances of 20 kpc behind the galaxy. 

Another feature that appears in all simulations is 
the flittering of the outer edge of the galaxy. The KH- 
instability pulls low density tongues along the upstream 
side of the remaining gas disk from the centre towards the 
edge. There these pieces of gas are exposed to the full wind 
and blown away (see e.g. last two panels in Fig. lll() . This 
procedure repeats continuously and leads to an oscillation 
in the radius of the gas disk rdisk(i)- 

The evolution of the density in the galactic plane and 
the surface density profile are shown in Fig. 1141 for a 
stronger and for a weak subsonic wind. As long as a sub- 
stantial part of the gas disk remains, neither the density 
profile nor the surface density in the inner part change 
much. However, we observe a compression of the outer 



layers on the upstream side of the gas disk (see Sect. 15. 5f) . 
For low ram pressures only the very outer layers are com- 
pressed, for stronger winds the compression goes deeper 
towards the galactic plane. 

5.3. Varying the wind parameters 

In this Sect, we show how the results for the gas disk ra- 
dius r(iisk(i)i the gas mass inside the original disk region 
Morig(i) and the bound gas mass Mbnd(0 depend on the 
wind strength. We also show in separate plots (in order to 
maintain clarity) how much bound gas can be found out- 
side the original disk region (that is Mbnd(0 ~ -Aforig(O) 
and how much fallen back gas is contained inside the orig- 
inal disk region (this is Mtb(t), see Sect. I4.5"T|) . We scan 
the wind parameter space for three different cases: 

(a) The massive galaxy in different ICM winds with ICM 
temperature Tjcmi (corresponding to an ICM sound 
speed of 1000 km s'^). 

(b) The massive galaxy in different ICM winds with ICM 
temperature Ticm2 (corresponding to an ICM sound 
speed of 562 kms~^). 

(c) The medium galaxy in different ICM winds with ICM 
temperature Tjcmi (corresponding to an ICM sound 
speed of 1000 km s^^). 

As a representative case, the results of case (a) are 
shown in Figs. 1151 [TCI and El Corresponding plots 
for cases (a) and (b) are very similar. For clarifica- 
tion we apply a time-averaging to most rdisk(i) curves, 
smoothing the oscillations that are due to the flitter- 
ing edge of the galaxy (see end of Sect. 15. 2|) . We av- 
erage rdisk(i) over 100 Myr for the cases with pram > 
1000 cm-3 km^ s-2 (p^am > 100 cm^^ km^ s"^ for the 
medium galaxy) and over 200 Myr for the lower ram pres- 
sure cases. The smoothing is applied neither to the ini- 
tial 250 Myr (500 Myr) for the stronger ram pressures 
(for the weaker ram pressures), nor to the curves for 
Pram = 10 cm~'^ km^ s~^ . Also in each group of constant 
ram pressure we leave one curve oscillating. The amplitude 
and the frequency of the oscillations are similar within a 
group. 

The first thing to notice is that the mass and radius of 
the remaining gas disk depend mainly on the ram pressure 
with higher ram pressures resulting in smaller disks. For 
the strongest ram pressures even the complete gas disk is 
stripped. At second glance one can see that for a given 
ram pressure the gas loss is stronger for subsonic winds 
than for supersonic ones. 

The phases of the stripping process as they were ex- 
plained in Sect. 15.21 are clearly visible in the rdisk(i), 
Mo^ig{t) and Mbnd(^) curves. The bending of the outer 
gas disk during the instantaneous stripping phase results 
in a first strong decrease of rdisk(0 ^^.^ MDrig(i)- We define 
the duration of this phase Tpush to last from the moment 
the flow reaches the galaxy to the first local minimum in 
fdis\^{t)- Fig. 1271 shows the dependence of Tpush on ram 
pressure, this point is discussed in Sect. 15.7.11 After the 
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time / Myr 

Fig. 15. Gas disk radius rdisk(i) (top panel), mass Morig(i) 
in the original disk region (middle panel) and bound gas 
mass Mbnd(i) (bottom panel) of the massive galaxy in 
ICM winds with Ticmi for different ram pressures Pram 
and Mach numbers. Curves for the same Pram group to- 
gether. Each group is labelled accordingly with its Pram in 
units of cm~'^km^s~^. The Mach number of each curve 
is coded by the line style, see key in the middle panel. 
The right axes are labelled in units of the initial radius or 
mass. 

In the top panel all curves except the ones for Mach 0.8 
are shown in the smoothed version (see text) for clarity. 
The frequency and amplitude in rdisk(i) are comparable 
for all runs with the same pram- 

The vertical bars show where we consider the intermediate 
phase to be finished. The two thin vertical lines mark the 
moment when the flow reaches the galaxy {to = 40 Myr) 
and to + ^switch = 60 Myr. 



instantaneous stripping phase, the intermediate dynamic 
phase follows with the falling back of some gas. This can 
be seen in the slight increase of Morig(i) just after the 
initial decrease, and of course in the plot showing M{h{t) 
directly fFig. I17|l . Although the gas disk always reacts to 
the wind in the initial instantaneous stripping phase, only 
for stronger winds this also results in unbinding a sub- 




200 400 600 800 1000 1200 1400 1600 1800 2000 
time / Myr 

Fig. 16. Mass of bound gas outside the original disk region 
as a function of time; for the massive galaxy in different 
ICM winds (Ticm)- There is one labelled panel for each 
ram pressure. Inside the panels the Mach numbers of the 
wind is coded by the line style (see key). The right axes 
are labelled in units of the initial radius or mass. 
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Fig. 17. Mass of fallen back gas Mfb(t)(= galactic gas 
that has left the original disk region and fell back to that 
region) as function of time; for the massive galaxy in dif- 
ferent ICM winds (Ticmi). There is one labelled panel for 
each ram pressure. Inside the panels the Mach numbers 
of the wind is coded by the line style (see key). The right 
axes are labelled in units of the initial radius or mass. 



stantial amount of gas from the galactic potential. The 
duration of the intermediate phase, this is until most of 
the bound gas outside the disk region has vanished, takes 
much longer for weak winds than for strong ones. For ram 
pressures of 10cm~'^km^ s^^ this phase is not finished at 
the end of the simulations. For all but the very weakest 
ram pressure, during the first few 100 Myr between 20 and 
40% of the original gas disk mass linger outside the orig- 
inal disk region but are still bound. The amount of gas 
that falls back to the original disk region is of the order 
of 5 to 10% of the original gas disk mass. 
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Fig. 19. Density contours at 70Myr for three different 
gas disk tliicknesses (6gas — 0.2, 0.4, 0.8 kpc). Gas in the 
shaded region is bound to the galactic potential. The wind 
was the same in all cases (pram — 1000cm~'^km^s~^, 
TicMi, Mach 0.8). The method of radius measurement is 
illustrated (see Sect. I4.5TT|) . 

For the weakest ram pressure (pram = 
10 cm~'^ km^ s~^) the galaxy loses only a few per- 
cent of its gas, although also in this case the radius of the 
gas disk is truncated. 

5.4. Influence of the wind initialisation 

We test how our results are influenced by the initialisa- 
tion of the simulations (see Sect. \A.'2]i . In addition to the 
standard case ^switch — 20 Myr we rerun a few represen- 
tative cases with tswitch — 70 Myr and ^switch = 150 Myr. 
The resulting mass and radius of the gas disk for two of 
these winds are shown in Fig. 1181 Again we smooth the 
curves for rdisk(0 (see Sect. l5.3jl . Different ^switch result in 
expected systematic differences during the initial phase. 
The decrease in mass and radius due to the instantaneous 
stripping is fastest for the shortest ^switch- After that rdisk, 
Aforig and Mbnd agree well despite the different tswitch ■ The 
cases with the longest ^switch tend to retain a little more 
mass during the intermediate phase. For weaker Pram the 
agreement between cases with different iswitch is not as 
tight as for strong Pram, which is due to the longer dy- 
namic (chaotic) intermediate phase for weaker Pram- The 
qualitative behaviour of the amount of bound mass outside 
the disk region and of the fallen back mass is independent 

of ^switch ■ 

5.5. Influence of the vertical structure of the gas disk 

In order to study the influence of the vertical structure 
of the gas disk we perform some simulation runs with ex- 
ponential and flared disks of varying scale heights in rep- 
resentative wind cases (see Table 13). Fig. 1191 shows snap- 
shots at the same moment during the instantaneous strip- 
ping phase for three exponential disks with different scale 
heights. Despite their different thickness, all three cases 
look very similar. The resulting radius and mass curves are 
compared in Fig.lJOlfor two of the winds. It does not make 
sense to use Morig for this comparison, as for gas disks with 
different thicknesses and shapes also the original disk re- 
gions differ strongly. Hence we use Mcyi here, the mass 



inside the cylinder with r < 26 kpc and \z — Zgai| < 5 kpc. 
The different starting points of the curves for the bound 
mass are due to the different disk thicknesses and external 
pressures (see Sect. 13. 2.'^ and Fig.EJ; for thicker disks the 
outer layers become unbound. The radius curves are again 
smoothed as explained in Sect. 15.31 

The Tdisk as well as Mcy\ and Mbnd are nearly indepen- 
dent of the vertical structure of the gas disk. For the higher 
ram pressure of 1000 cm^"^ km^ s^^ the results for the dif- 
ferent gas disks agree excellently. For the lower ram pres- 
sure (100 cm~'^ km^ s~^) the scatter is a bit larger as the 
galaxies stay longer in the dynamic intermediate phase. 
The difference in the radius curves during the initial in- 
stantaneous stripping phase are due to the method of mea- 
surement, as is demonstrated in Fig. El We conclude from 
this test that the success of RPS is independent of the ver- 
tical structure of the gas disk. 

We compare the compression of the upstream side of 
the gas disk for different scale heights and winds by show- 
ing the evolution of the vertical density profile through 
the gas disk at r = 0.23 kpc in Fig. [2] The compression 
of the upstream side of the gas disk results in a steep- 
ening of the density profile of the outer layers. In the 
compressed region the density profile can be fitted with 
an exponential function of a smaller scale height than the 
original bgi^. The scale height in the compressed layers is 
about 0.12 kpc independently of the wind and the origi- 
nal fogas, but for stronger ram pressures the compressed 
layer reaches deeper towards the galactic plane. The de- 
pendence of the scale height of the compressed layer on 
resolution is demonstrated in the bottom line in Fig. 1211 
it is only a weak dependence. This plot shows vertical pro- 
files at t = 100 Myr and in the time intervals 300-400 Myr 
and 700-800 Myr. We chose to show time intervals to allow 
averaging over the intrinsic oscillations in the profiles. In 
any way, if the scale height of the compressed layer was set 
by the resolution, it should vary from thick to thin disks. 
For the medium galaxy the profile of the outer layers is 
steepened to a scale height of about 0.09kpc. 

5.6. Cross-comparisons between runs with different 

TlCM 

The results from Sect. 15.31 suggest that the effect of the 
stripping does not depend on the ram pressure alone 
but also slightly on the Mach number. To investigate 
this point further we compare runs with different ICM 
temperatures. E.g. for Tjcmi the velocity of 800knis~^ 
is equivalent to Mach number 0.8, whereas for Ticm2 
the same velocity corresponds to Mach number 1.423. 
The velocity that belongs to Mach number 0.8 for Ticm2 
is 450 km s^^. If indeed the Mach number is the deter- 
mining parameter, then for this example the result of 
the run with (Tjcmi, 800 kms^^, Mach number 0.8) 
should be the same as the result of the case 
(TicM2, 450 kms^^, Mach number 0.8) rather than 
(TicM2, 800 km s^S Mach number 1.423). We show this 



16 



Elke Roediger and Gerhard Hensler: Ram pressure stripping of disk galaxies 




200 400 600 800 1000 1200 1400 1600 1800 2000 200 400 600 800 1000 1200 1400 1600 1800 2000 



time / Myr time / Myr 

Fig. 18. Comparison of stripping radius rdiskit), niass in original disk region Aforig(i) and bound mass Afbnd(0 f'^'" 
different ^switch, for two representative winds. We marked the moment where a sound wave starting at the inflow 
boundary reached galaxy (to), and to + tswitch for the three different Switch by the four thin vertical lines. For further 
explanations see Fig. ^1 



comparison for rdisk(i), Moiig{t) and Mbnd(i) in Fig. [221 
for some cases. In the middle panel, where the velocity 
is subsonic for TicMi and supersonic for Ticm2, indeed 
runs with the same Mach number rather than with 
the same velocity correspond to each other. The same 
seems to hold for the left panels, in the pure subsonic 
range, although not as clear as in the transition from 
the subsonic to the supersonic regime. In the supersonic 
range (right panels) no systematic difference is seen. As 
qualitative issues we conclude that (i) for a given wind 
density and velocity, RPS is slightly more effective if this 
wind is subsonic (compared to supersonic winds); and 
(ii) that in the subsonic regime winds with smaller Mach 
numbers are a bit more effective than winds with higher 
Mach numbers. Nevertheless, the Mach number is only 
a secondary parameter, the main parameter is the ram 
pressure. 



Moreover, we do not want to over-interpret the dif- 
ferences between the subsonic and supersonic cases, as 
we used a single-fluid description. In contrast, the ICM 
is highly ionised and actually consists of a proton and 
an electron fluid. With respect to the proton fluid, the 
galaxy's velocities are transonic, but with respect to the 
elec tron fluid they still move subsonically. ^Po rtnoy et alJ 
()l99 3) studied the difference between the single-fluid and 
a two-fluid description for a spherical galaxy (with gas re- 
plenishment) and found no difference for the amount of 
gas remaining inside the galaxy. 



5.7. Comparison of analytical estimate and numerical 
result 

5.7.1. Instantaneous stripping 

We have found that during the instantaneous stripping 
phase the outer part of the gas disk is displaced, but not 
immediately unbound. It takes a while until the displaced 
gas is also unbound. As the analytical estimate in Eq. ^ 
just compares the forces working on the gas disk, it tells 
how much gas will be unbound, but it cannot quantify 
how long this process may take. As the analytical estimate 
states how much gas will be unbound from the galaxy, 
we also need to extract this quantity from the simula- 
tions. Hence we measure the radius and the mass of the 
remaining gas disk at the end of the dynamic intermediate 
phase, i.e. when no substantial amount of bound gas out- 
side the disk region is left. We have marked this moment 
in Fig. 1151 by vertical bars. For the lowest ram pressure 
(Pram = 10 cm~^ km^ s~^ ) the intermediate phase is not 
completed during the simulation run time, here we use 
the final values for the comparison. 

As we pointed out in Sect. I5.im there are several ver- 
sions how the analytical estimate can be done, therefore 
we discuss the following recipes: 

(a) Use Eq. ^1 to compute the restoring force while using 
the maximal potential gradient inside the original disk 
region (marked by "E/in disk" in Figs. 1231 and I24() . 

(b) Use Eq.^lwith Az — 150 pc and use the maximal po- 
tential gradient inside the original disk region (marked 
by V/in disk" in Figs. [23 and EH). 
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Fig. 20. Comparison of stripping radius r^iskit), mass in 
different vertical structures, for two representative winds, 
smoothed as explained in Sect. 15.31 
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original disk region A/orig(t) and bound mass Mbnd(0 f'^'" 
For further explanations see Fig. 1151 All radius curves are 



(c) Use Eq. ^] with Az = 150 pc and use the maximal 
potential gradient found for any z (marked by "p/ all 
z" in Figs. 1^ and Eljl . 

(d) Compare the thermal pressure in the galactic plane 
and the ram pressure (Eq. 1161 marked by "pressures" 
in Figs. 113 and EH) 

Version (c) is between (a) and (b) ; for higher ram pressures 
it is indistinguishable from version (b), as there smaller 
radii are concerned where the steepest potential gradient 
is found inside the disk region. Also the result from the 
pressure comparison (d) has a very similar shape to (a) 
and (b) and is between those two. As already mentioned 
in Sect. l5.lTTl the estimates according to versions (a) and 
(b) are identical despite a shift in horizontal direction for 
exponential disks. These two versions represent the two 
extremes in the uncertainty how the gravitational restor- 
ing force /grav shall be computed - either simplifying to a 
thin disk (a) or computing /grav for a single cell only (b). 
This uncertainty can be parametrised by a factor of a in 
the calculation of /grav , modifying Eq. ^1 to 



I /grav (2, r) I 



dz 



(z,r)I]gas(r) 



(27) 



Choosing a = 1 corresponds to version (a), choosing a = 
0.16 recovers version (b), but a can be adjusted to fit the 
numerical results. 

Before we do any fitting, we show the pure analytical 
estimates (versions (a) to (d)) and the numerical results 
in Fig. 1231 for the massive galaxy and in Fig. for the 
medium galaxy. The numerical results are indeed in the 
range suggested by the analytical estimates. 



We also want to consider the Mach number depen- 
dence of the numerical results again. So far we calculated 
the ram pressure according to Prum — "-icmw^cm fo^ 
supersonic cases as well. However, the velocity and den- 
sity behind the formed bow shock differ from the values 
on its upstream side. It may be more appropriate to cal- 
culate the ram pressure for supersonic winds with the 
values of p and v at the downstream side of a straight 
shock with the corresponding Mach number. The com- 
ponent of the momentum density picuViCM perpendicu- 
lar to the shock front is conserved across such a shock, 
whereas the velocity at the downstream side is reduced 
(and the density increased) by a certain factor accord- 
ing to the Rankine-Hugoniot conditions. Then also the 
true ram pressure would be diminished by the same fac- 
tor as V. For Mach numbers of 1.423 and 2.53 (at the 
upstream side) the velocity at the downstream side is re- 
duced by factors of 0.62 and 0.36, respectively. In Figs. 1251 
and 1261 we replot the numerical results with this correc- 
tion, that is shifting the points for the supersonic results 
down to the corrected ram pressures. In addition, we re- 
peat the analytical estimate from the pressure compar- 
ison (version (d), marked with "pressures" in Figs. E5l 
and I26|l and plot the estimate from the force compari- 
son with an adjusted a (see Eq. 1271 marked with "forces, 
a = 0.25" in Figs. ES] and EHI)- For a = 0.25 the numerical 
result and the analytical estimate agree well. In contrast 
to the uncertainty expressed by a in the force comparison, 
there is no uncertainty in the pressure comparison. And 
indeed the (Rankine-Hugoniot-corrected) numerical result 
and the estimate from the pressure comparison are in good 
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Fig. 22. Comparison of stripping radius rdisk(i), mass in original disk region Morig(t) and bound mass Mbnd(i) for 
selected runs with pram — 1000 cm~'^ km^ s" 
explanations see Fig. El See also Sect. 15.61 



(upper panel) and PramlOO = cm km s (lower panel). For further 



agreement. With a = 0.25, the estimate from the force 
comparison was shifted to match the result from the pres- 
sure comparison. The point for the lowest ram pressure 
and Mach 2.53 in Fig.l^Olthat lies above the analytical esti- 
mate is explained easily because the dynamic intermediate 
phase for this run is not finished at the end of the simula- 



tion. We note that the correction of the ram pressures for 
the supersonic cases according to the Rankine-Hugoniot 
conditions cannot describe the situation completely, be- 
cause it would mean that also in the supersonic range the 
results should depend slightly on the Mach number. This 
is not found in the simulations (see Sect. 15. 6f) . 
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Fig. 21. Evolution of vertical profiles a,t r — 0.23 kpc for 
different -winds and gas disk thickness. We used an ex- 
ponential disk for all cases, the scale heights -were 6gas = 
0.2 kpc (thin), &gas = 0.4 kpc (standard) and &gas = 0.8 kpc 
(thick). We fitted an exponential function to the steep- 
ened part of the profile, the scale height for this fit -was 
approximately 0.12 kpc for all cases. In the bottom line 
the variation of the vertical profiles -with resolution for 
three different time steps/time intervals is demonstrated 
(colour version in online paper). 



We can also compare the predicted duration for the in- 
stantaneous stripping -with the numerical result. The pre- 
diction dra-wn from Eq. 1221 refers only to the replacement 
of the outer gas disk. The numerical result for the replace- 
ment is given by the time of the first local minimum in the 
'"disk curves in Fig. 1151 The trends in the numerical result 
and the analytical estimate sho-wn in Fig. 1271 agree -well, 
also the absolute values agree roughly. We also sho-w the 
duration of the intermediate phase, that is the time until 
all gas that is not bound tightly enough is truly unbound. 
We find that this time is approximately 10 times as long 
as the instantaneous stripping phase. 
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Fig. 23. Comparison of analytical estimate and numer- 
ical result for the stripping radius and the remaining 
gas mass as a function of ram pressure, for the massive 
galaxy. We sho-w four versions of the analytical estimate, 
see Sect. l5.7Jl The Mach numbers for the numerical result 
are coded -with different symbol shapes. Open symbols are 
for cases with Tjcmi, solid symbols of the same shape for 
cases -with Ticm2- 
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Fig. 24. Same as Fig. |231 but for the medium galaxy. 



5.7.2. Continuous stripping 

The estimate of the mass loss rate during the continu- 
ous stripping phase (see Sect. 15.1.21) is only a rough one. 
The mass-loss rates are expected to be of the order of one 
Mq yr^^. Fig. 1281 demonstrates the mass-loss rates de- 
rived from the simulations by fitting a linear function to 
Mhnd{t) during the continuous stripping phase. Because 
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Fig. 25. Comparison between the analytical estimates and 
the numerical result for the massive galaxy with Rankine- 
Hugoniot-corrected ram pressures for the supersonic cases 
(see text). Open and solid symbols are for cases with 
TicMi and Ticm2, respectively. The symbol shape codes 
the Mach number. The analytical estimate for the force 
comparison (solid line) is fitted to the numerical results 
via the parameter a (see also text). The estimate from 
the pressure comparison is repeated unmodified. 
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Fig. 27. Duration of the instantaneous stripping phase 
and intermediate phase as a function of ram pressure. 
Large symbols are for the massive galaxy, small symbols 
for medium galaxy. Also the analytical estimate for the 
instantaneous phase according to Eq. 1221 is shown, the 
solid line is for the massive galaxy, the dashed one for 
the medium galaxy. 



03 



m 

(0 

o 

(0 

in 
ca 
E 




-'ram="ICM ^ICM 



in cm 



Fig. 26. Same as Fig. but for the medium galaxy. 



Fig. 28. Mass loss rates of the continuous stripping phase 
as a function of ram pressure. The Mach numbers are 
coded by different symbols. Symbols of group (I) (large 
symbols) are for the massive galaxy in winds with Tjcm i , 
group (II) (medium-size symbols) for the massive galaxy 
in winds with Ticm2i and group (III) (small symbols) for 
the medium galaxy in winds with Ticmi- We also show 
the rough analytical estimate (see Eq. 12^1 in Sect. 15. iT^ . 



est mass-loss rates, we cannot find any further trends. The 
analytical estimate could give a rough impression only. 



the runs with Pram — 10 cm^^ km^ s~^ do not reach the 
continuous phase during the runtime, for those runs we fit 
the time range between 1 and 2 Gyr. From the simulations 
we can derive mass-loss rates between 0.1 and 2 Mq yr^^. 
Except that the highest ram pressures produced the high- 



6. Discussion and conclusions 

6.1. Influence of the ICM pressure 

What happens to the gas disk of a spiral galaxy in 
a high pressure external medium? From the compari- 
son of typical pressures in the disks of isolated galax- 
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ies (10~^^ to 10~^^erg cm^^) and typical ICM pressures 
(e.g. 1.4 • lO^^'^ erg cm^'^ for tiicm = and 
TiGM = 10^ K) alone it is obvious that normal disk galax- 
ies cannot exist in cluster environments, no matter if they 
move or not. If for an isolated galaxy the external pres- 
sure is increased gradually, the disk has two alternatives 
to react to that. Either, if the density distribution should 
not change much, it has to heat up in order to reach 
higher pressures so that it can maintain pressure equilib- 
rium with the surrounding gas. This is what we assumed 
implicitly with our initial conditions. Or, if cooling pro- 
hibits a substantial increase of the temperature, the disk 
will be compressed. However, this is only the crude ba- 
sic idea. The fate of a more realistic multi-phase gas disk 
including effects like star formation and feedback may be 
much more complicated. Due to the compression the for- 
mation of molecular clouds and stars may be enhanced, 
which in turn could increase the amount of heating by an 
increased supernova rate and stellar winds. The observa- 
tional fact that even Hl-deficient cluster galaxies contain a 
rather normal am ount of molecular gas l|Kennev fc Yound 
ll989HBosem et al...l997;.Cason et al.lll998|) and show nor- 
mal to enhanced star formation rates in their inner part 
ijKoopmann fc Kennevlll998l l2nn4alhD supports this pic- 
ture. It would be worth to study the effect of the external 
pressure on the gas disks of spiral galaxies in a separate 
work. 



In how far does our setup of the initial model (fixing 
the density distribution of the gas disk and calculating the 
temperature distribution from hydrostatic equilibrium), 
which leads to high disk temperatures, bias the results 
derived in this paper? The concern is that due to the high 
temperature (and hence higher thermal energy) our gas 
disks are not bound as strongly as cooler disks in isolated 
galaxies and hence are stripped easier. In this case we 
would overestimate the stripping efficiency. The alterna- 
tive way of initial setup would be to fix the temperature 
inside the gas disk (e.g. determined by the heating-cooling 
balance) as well as its mass and calculate the density dis- 
tribution from the hydrostatic equation. The resulting gas 
disks would become thinner and have higher densities as 
the external ICM pressure increases. In varying the gas 
disk thickness and shape (see Sect. 15.511 we have already 
checked if such thin disks are harder to strip by ram pres- 
sure, but we did not find them to be more resistant to 
RPS. We also performed test runs that started with a 
galaxy resembling an isolated one (concerning disk pres- 
sure and temperature), but put into a high-pressure ICM. 
As this galaxy was not in hydrostatic equilibrium with the 
ICM, it was compressed immediately and then stripped 
like its equilibrium analogue. Hence we conclude that our 
results are not biased strongly by the initial pressure and 
temperature distribution. 



6.2. Morphological features of stripped galaxies 

In agreement with previous work, our simulations sug- 
gest that disk galaxies suffering (approximately) face-on 
RPS should have truncated disks, where the radius of the 
disk is set by the ram pressure and the galactic poten- 
tial. Moreover, if they are in the instantaneous stripping 
phase or the dynamic intermediate phase, stripped gas 
could be found at distances as far as 20kpc behind the 
disk. How long such structures exist depends at least on 
Pram- In how far this stripped gas actually remains visi- 
ble, is dispersed or cools, forms stars or is affected by heat 
conduction, our simulations cannot predict as they do not 
include such processes. Due to the homogeneous gas disk 
used in our simulations, filaments emerging from the disk 
are found mainly at the disk e dge. Such examples ar e in- 
deed observed (e.g. NGC 4522 :IVollmer et al. ( 2004b) and 
references therein, NGC 4402: ICavatte et al.i (4990.1199^ : 
ICrowl et al.l i2no4\) The filaments found in NGC 4569 
llTschokc ct al. 20 0lUVollmer et al.ll2004at iKennev et all 
I2OO4 .Romans et al.ll20o'^ seem not to emerge from the 
edge but more from the centre of the disk. This may be 
due to activity in the galactic nucleus or due to RPS of 
an inhomogeneous gas disk. 

6.3. Mixing with the ICM 

We tried to infer to what extend ICM gets mixed into 
the remaining gas disk. Such mixing would change the 
angular momentum and metalicity of the gas disk. To in- 
vestigate this point we used two more recipes to compute 
the mass of the gas disk or the bound mass. In contrast to 
summing over the "coloured" gas (the gas that has been 
inside the galactic disk initially) in the disk region, we 
also summed up all gas in the disk region with a temper- 
ature below 5 • 10^ K (about the virial temperature of the 
massive galaxy). This prevents that we include the ICM 
that is streaming freely through the already swept disk re- 
gion, but it would include all ICM that has already mixed 
with the gas disk. On the other hand, we sum over all gas 
in the region of the bound gas instead of summing only 
over the "coloured" gas there. This includes all ICM that 
is situated in this regions. We find, however, no differ- 
ences between the versions that should include or exclude 
the ICM, so we conclude that no ICM is accreted by the 
galaxy. 

6.4. Time dependent winds 

During its passage through a cluster a galaxy does not 
experience a constant wind. Galaxies on more radial or- 
bits fly through the high density cluster core. With a 
typical velocity of 1000 kms^^ a galaxy needs 1 Gyr to 
cross the inner Mpc of the cluster. lyoTlmer et alJ l)2001bl) 
calculated the time dependent ram pressure for a galaxy 
falling into the Virgo cluster (their Fig. 3) and found that 
typically the ram pressure exceeds 1000 cm^"^ km^ s~^ for 
about 100 Myr before and after the core passage, and still 
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exceeds 300 cm~^ km^ s~^ about 250 Myr before and after 
the core passage. According to our simulations these times 
are definitely long enough to truncate the radius of the gas 
disks. However, the duration of the high ram pressure may 
be too short to unbind all this gas. Even in our simulations 
with constant winds several lO^AfQ fall back to the disk 
region, and for quite a while about IO^Mq of gas that is 
still bound lingers behind the galaxy. If the ram pressure 
abates too early after the core passage all this gas should 
fall back to the disk. From our simulations we cannot in- 
fer in how far and how fast the gas disk resettles after the 
core passage. 

6.5. Galaxies in cluster outskirts and galaxy groups 

In our simulations we found that even mild ram pres- 
sures of about 100 cm^'^ km^ s~^ and less already trun- 
cate the gas disk of massive galaxies to about 15kpc. 
Even if low ram pressures do not unbind much gas, 
they still displace the outer parts of the gas disk. The 
wind itself and back-falling gas disturb the disk, send 
shock wav es through it which could trigger star forma- 
tion (e.g. iFuiita fc~Nagashinialll999f l. We therefore ex- 
pect that even at larger distances from the cluster centre 
the gas disks of most galaxies should be "pre-processed" 
(see also lFuiitalEo041 . This pre-processing of galaxies in 
cluster outskirts and groups could explain the HI deficient 
galaxies at large dis tances to the cluster centre found by 
ISolanes et all l)200j) . RPS in low density environments 
may be more important than thought so far. 

6.6. Inclination 

We have studied only cases where gala xies move face- 
on. W hat can we learn for inclined cases? iMarcolini et alJ 
l(200 3'l found that the inclination angle matters only for 
medium ram pressures. This can be understood also in the 
context of our simulations. For very weak winds hardly 
any gas is unbound from the galaxy by the pushing effect, 
the main mechanism is the continuous stripping. For this 
mechanism the inclination angle is not important, it works 
on the whole surface. For very strong winds the complete 
disk is lost independently of inclination. Only for the in- 
termediate winds the inclination decides how much mass 
can be pushed out in the beginning, hence for heavily in- 
clined cases less gas is lost in this phase. After that the 
continuous stripping continues at the stage where the pre- 
vious phase has ended. Also lOuilis et al.l l)2000|) agreed 
that only strict edge-on galaxies can retain more gas, in- 
termediately inclined galaxies are stripped similarly to 
fa ce-on ones. Th i s is als o in agreement with the results 
of IVollmer et al.l l(2001hl) . who constructed a generalised 
formula relating the remaining gas mass with the ram 
pressure and the inclination angle. However, their sticky- 
particle s imulations canno t mod el the turbulent viscous 
stripping. ISchulz fc Striickl l|200lh found that for inclined 
galaxies less gas lingers in the lee of the disk, hence our 



results concerning back-falling gas may change in such sit- 
uations. However, the unbinding of the gas after it has left 
the disk region should take some time also for inclined 
cases. 

In addition to the restriction to the face-on case, a 
2D code is restricted to rotationally symmetrical cases. 
However, real systems are not symmetrical, and an asym- 
metry may change the behaviour of turbulence and insta- 
bilities. 

6.7. Summary 

We performed a large set of numerical simulations of the 
face-on RPS of disk galaxies in constant ICM winds. We 
find that even massive galaxies can lose a substantial part 
of their gas disks, and can be stripped completely in clus- 
ter centres. Our simulations also covered low density en- 
vironments like cluster outskirts, and we find that there 
galaxies are truncated to 15 to 20 kpc. Back-falling gas and 
the compression from the flow could trigger star forma- 
tion, hence we expect galaxies to be pre-processed in clus- 
ter outskirts and groups as well. Furthermore, the models 
show that the stripping proceeds in three phases: 

— Instantaneous stripping phase: the outer disk is pushed 
out of its initial position (bent towards the downstream 
side), but it is only partly unbound. The main radius 
reduction happens in this phase. This phase is short 
(20 to 200 Myr), its duration as a function of ram pres- 
sure can be read from Fig. 1271 

— Intermediate phase: rather dynamic, the gas outside 
the disk region partially falls back to the disk (about 
10% of the original gas mass), the rest is completely 
stripped. This may repeat in cycles. This phase lasts 
about 10 times as long as the instantaneous stripping 
phase, it is finished when no substantial amount of 
bound gas is left outside the disk region. The radius 
does not change much during this phase, but the bound 
gas mass and also the gas mass inside the disk region 
do. The mass and radius of the remaining gas disk at 
the end of this phase can be predicted best by com- 
paring the thermal pressure in the disk plane and the 
ram pressure (for details see Sect. 15. I'll Eg. I16|l. The 
analytical estimate and the numerical results for the 
stripping radius and the remaining gas mass are sum- 
marised in Figs. 1211 and [SB] for the massive and the 
medium mass galaxy, respectively. 

— Quasi-stable continuous stripping: constant small mass 
loss rate of about IMq yr^^, slow decrease of the disk 
radius. 

The main parameter that sets the mass and radius of the 
remaining gas disk is the ram pressure, but there is a slight 
dependence on Mach number in the sense that in the sub- 
sonic regime and at the transition from the subsonic to 
the supersonic regime ICM flows of lower Mach numbers 
result in a bit stronger gas loss than flows with higher 
Mach numbers. In how far this result holds for a two-fluid 
treatment of the ICM is unclear. 
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Pram 


Mach number 


viscosity parameter 


(cm"^ km^ s~^) 




C2 (see Sect. HI 


1000 


0.8 


1,2,6 




2.53 





Table A.l. List of simulations to study the influence of 
the artificial viscosity. Common parameters for all runs 
are: massive galaxy with standard exponential disk, ICM 
wind with Ticmi (cicm — 1000 km s^^). 



The mass loss and shrinking of the radius depend only 
on the surface density of the gas disk, but not on its ver- 
tical structure (thickness, flared or exponential). 

We observe a compression at the upstream side of the 
remaining disk, but we do not find an increased surface 
density. 
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Appendix A: Influence of artificial viscosity 

In our simulations we observed the long phase of turbu- 
lent viscous stripping. The code uses an artificial viscosity 
mainly to be able to handle shocks. A concern may be 
that due to this artificial viscosity the code cannot handle 
the KH-instability correctly and hence biases the results of 
the viscous stripping phase. We performed test runs (see 
Table IX!T)l to investigate the infiuence of the artificial vis- 
cosity by varying the viscosity parameter C2 (see Sect.EJ. 
For the results for the disk radius, mass and bound mass 
the differences between runs with various viscosities are 
negligible. We show one case in Fig. lA.ll Fig. IA.2I com- 
pares snapshots at t = 200 Myr for the test cases with 
varying viscosity. As expected, for the high viscosity pa- 
rameter, the surface instabilities are softer, whereas for in 
the low viscosity case (C2 = 1) the features are sharper. 
However, almost all features can be identified for all vis- 
cosities. We conclude that the use of the artificial viscosity 
does not bias our results. 

Appendix B: Influence of resolution 

In order to test the influence of the resolution we have 
performed test runs with resolutions 50 pc, 100 pc, 150 pc 
and 200 pc. 

Fig. IB. II compares the radius rdisk and mass Mdisk 
of the remaining gas disk for the standard run (mas- 
sive galaxy, exponential gas disk, 6gas = 0.4 kpc, ticMi, 
Pram = 1000 cm""^ km^ s~^, Mach number 0.8) for the four 
different resolutions. Also the the mass of bound gas Afbnd 
and the mass of fallen-back gas Mfb are shown. The differ- 
ences for all quantities are negligible. The small differences 
in Mbnd arise because in cases with higher resolution the 
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Fig. A.l. Influence of the viscosity on the evolution of the 
gas disk radius, mass and fallen-back mass. The value of 
the viscosity parameter C2 (see Sect.El) is indicated in the 
key. For further parameters of these runs see Table EH 
This plot is for the runs with Mach number 0.8. 
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Fig. A. 2. Comparison of the density distribution (at t = 
200 Myr) for different viscosities. The plot style is the same 
as in Fig.^] The top panels are for Mach number 0.8, the 
bottom panels for Mach number 2.53. For further param- 
eters see Table IXm Each panel is labelled with the appro- 
priate viscosity parameter C2 (see Sect.|2). The rectangle 
corresponds to 10 x 10 cells. 



the stripped gas can be compressed into smaller volumes, 
reaching higher densities locally. Such dense clouds are 
more likely to be bound than the clouds with lower den- 
sities. 

In Fig. IB. 21 we show snapshots a,t t — 200 Myr for the 
four different resolutions. The strongest difference occurs 
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Fig. B.l. Influence of grid size and resolution on evolution 
of gas disk radius, mass and fallen-back mass. See text for 
parameters. 
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t^-Zgai) / Kpc 

Fig. B.2. Influence of resolution: Snapshots at < = 
200 Myr for four different resolutions as given in the plot, 
see text for parameters. The plot style is the same as in 
Fig. ^] The size of 10 x 10 cells is indicated by the rect- 
angle in each panel. 



for the stripped material. The coarser grids cannot resolve 
the small-scale fragmentation of the stripped gas evident 
in the highest resolution. However, we did not attempt 
to interprete the fate of the stripped material shown in 
the simulations, as we neglected e.g. thermal conduction. 



which might be play an important role in the evolution of 
the stripped gas. 
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